Interaction-assisted propagation of Coulomb-correlated electron-hole pairs in 

disordered semiconductors 
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A two-band model of a disordered semiconductor is used 
to analyze dynamical interaction induced weakening of local- 
ization in a system that is accessible to experimental veri- 
fication. The results show a dependence on the sign of the 
two-particle interaction and on the optical excitation energy 
of the Coulomb-correlated electron-hole pair. 

PACS numbers: 73.2Q.Dx, 72.15.Rn, 78.66.-w 



The problem of two interacting particles (TIP) in a 
random potential is an excellent paradigm for the general 
question of the interplay of disorder and interactions in 
many-body systems. First addressed in a 1990 paper by 
Dorokhov 0], the subject has been especially well studied 
since Shepelyansky's 1994 publication ||]. Considering 
the TIP localization length I2 , most authors obtain 
an interaction-induced increase I2 > h over the single- 
particle localization length li independent of the sign of 
the interaction, with 12/h ~ 11 and a = 1 or a = 0.65 
Here, h and I2 are measured in units of the lattice 
constant of a one-dimensional Anderson chain. Similar 
results have been obtained for TIP in a quasiperiodic 
chain Q . The independence of the predicted effect on the 
sign of the interaction is an especially intriguing feature. 

Early work approached this problem using a wide vari- 
ety of theoretical techniques and focused on establishing 
the existence of the TIP-effect, while more recent work 
has dealt with quantitative details like scaling behaviour 
and the influence of interaction range, strength and sign 
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Existing works comprise purely theoretical case stud- 
ies since the model of just two particles in a single band 
does not correspond to any real physical situation. Ob- 
viously, experimental study is needed to promote fur- 
ther understanding of the TIP problem, and would put 
the presently rather academic discussion on a firm phys- 
ical basis. Exploiting the fact that the coherent spatio- 
temporal dynamics of Coulomb-correlated electron-hole 
pairs is strongly influenced by the two-particle interac- 
tion, we show in the present paper that the TIP local- 
ization properties should be accessible to modern ultra- 
fast optical techniques. The corresponding spatial and 
temporal scales are on the order of sub-/im and 1 ps. 
Our numerical studies are based on integration of the 



Semiconductor Bloch Equations and include no a pri- 
ori assumptions about energy hierarchies or interaction 
matrix elements. The calculations were performed for a 
disordered ID semiconductor quantum wire, where local- 
ization effects are most important. Despite its simplicity, 
this model system contains already all essential ingredi- 
ents to describe the dynamics following optical excita- 
tion even in disordered systems of higher dimension [p^ . 
The model parameters have been given values that re- 
semble those of realistic disordered semiconductor quan- 
tum wires. Additional physical parameters (excitation 
energy, spectral pulse width, screening length, different 
masses of the two particles) allow the study of a wide 
variety of observable phenomena. 

In our numerical calculations we investigate the 
spreading of an electron-hole wave packet after local ex- 
citation by an optical pulse. Here the interaction is given 
by the long-range Coulomb potential which, besides pro- 
ducing bound states (excitons) near the edges of the ex- 
citation spectrum, also correlates the electrons and holes 
in the pair continuum. Previous theoretical studies of the 
spatial-temporal dynamics of wave packets formed from 
excitons show that their motion is rather limited in the 
presence of scattering |l^. This result is recovered by 
our present calculations. Here we focus our interest on 
the dynamics of optically generated wave packets in the 
pair continuum. We find that the excitation conditions in 
the presence of particle-particle interaction influence the 
carrier dynamics dramatically. In addition, and in con- 
trast to some previous claims in the literature, we find 
that the sign of the interaction has a pronounced effect 
on the spatio-temporal dynamics. 

We consider a ID array of sites i with diagonal dis- 
order in both the valence band {vb) and the conduction 
band (cb) . The site energies £vi and Ed corresponding to 
the vb and the cb, respectively, are randomly distributed 
over the interval [—W/2, W/2] and are uncorrelated. The 
nearest neighbor cb levels are coupled by the tunneling 
term — J^, the vb levels by — J". We use the Coulomb 
interaction in its monopole-monopolc form |l4[ with ma- 
trix elements 



U 



(1) 



which has been regularized in order to cope with the 
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pathological singularity in one dimension. The constant 
a has been chosen to be 5 times the lattice constant . 
The total Hamiltonian is written as H = Hq + Hj + He, 
where Hq describes the vb and cb band structure, Hj 
represents the semiclassical interaction with the electro- 
magnetic field Ei{t) in dipole approximation, and Hq 
defines the electron-electron interaction term ■ 

In the following, we assume a local initial excitation 
at the central site i = 0, which is modeled by setting 
fJ-i — ^i,o for the local dipole matrix element in Hj. The 
optical polarization pij (t) is obtained from the equation 
of motion for the polarization operator = diCj, which 
is coupled to the equation of motion of the electron and 
hole intraband quantities nfj = cfcj and hfj = dfdj, 

respectively, where the operators c^,Ci {df,di) describe 
the electron (hole) creation and annihilation operators 
at site i. The equation of motion for the expectation 
values Pij{t) and n^j^{t) is treated using the well-known 

Semiconductor Bloch Equations for pij {t) and n^j'' (t) in 
the real-space representation Detailed derivations 

of the Semiconductor Bloch Equations are given in p!^ ] 
and the textbook As we are interested in small 

excitation densities, we write only the equation for pij in 
the lowest (linear-response) order in the exciting field 

N 

dt Pt3 = - i ■ - 4 - V^j^p^j +i^{J'' Vti + J'' Pij) 

1=1 

+ i^i,E,it)S^j. (2) 

Using the conservation laws nfj = Y^jPij P;* and n''^ = 
TliiPjiPa valid in this lowest order |l3], we obtain the 
intraband quantities. 

Instead of studying a rather academic localization 
length which describes only the asymptotic behavior of 
wave functions, we calculate the experimentally more rel- 
evant participation number A(i) — (X^i^-fi)"^- Here na 
stands for either or n^^. With the packet localized at 
site 0, riii = 5iQ and A = 1, while for an excitation uni- 
formly extended over the sample of N sites, na = 1/N 
and A = A^. Our calculations were performed for chains 
containing N = 240 sites. Boundary effects can easily be 
identified in the temporal evolution of A and do not play 
any role as long as A < A^/2. All the data presented are 
free of finite-size effects. 

The transform-limited optical pulse is defined by its 
mean energy huj and the temporal width r of the gaus- 
sian envelope ~ exp{ — {t/T)^}. We define an excitation 
energy E^xc refered to the bottom of the (ordered) ab- 
sorption band, i.e. E^^c — hu — Egap. All results are 
given for r = 100 fs, which corresponds to an energetic 
width (FWHM) of 22 meV. 

To make contact with the previous work where two 
particles in a single band were placed initially at a single 
site, we first consider the situation of a symmetric band 
structure with J"^ = — j'' = 20 meV. The absorption 



spectra with and without Coulomb interaction are shown 
in Fig. |l| for the ordered case. The peak structure near 
the absorption is due to the excitonic resonances. Upon 
changing the sign of the Coulomb interaction, the bound 
state resonances are shifted from the bottom to the top 
of the absorption spectrum. 

As the dynamics of electrons and holes are the same for 
the assumed symmetric band structure, we restrict our 
discussion to the electrons. We first discuss the situation 
in the absence of Coulomb interaction. Fig. |^ shows the 
corresponding Ae{t) for two different disorder parame- 
ters W after excitation by a pulse at E^xc — 80 meV. 
The excitation is centered in the absorption spectrum as 
indicated in Fig. |l]. Ag (t) evolves exponentially with rise 
time less than 1 ps. Here and below, we take the satu- 
ration value as a measure of localization. As expected, 
it decreases rapidly with increasing disorder. We find 
A^M^~^-^asW^is varied over the range 40 meV to 240 
meV for J = 20 meV. A discussion of related exponents 
can be found in ||20[| . 

Fig. H contrasts the interacting and nonintcracting be- 
havior for two values of disorder and reveals three re- 
markable features, i) The interaction clearly leads to a 
reduction of the localization of the particles. We have 
carefully checked that the saturation value of Ae(i) at 
long times is not due to a finite size effect; values < A^/2 
are fully converged with respect to the sample size, ii) 
While the participation number in the nonintcracting sit- 
uation evolves exponentially and saturates quickly (< 1 
ps), the interacting wave packets evolve diffusively and 
reach their saturation values at much longer times, iii) 
The sign of the Coulomb interaction {U = ±1) has vir- 
tually no influence on the propagation of the particles in 
this case. The same is true if we apply a very short exci- 
tation pulse which spectrally covers the whole band. The 
spectral position of the central pulse frequency within the 
band is then completely irrelevant. In this situation the 
excited particle pair-wave packet is initially situated ex- 
clusively at site i = 0. These observations are not new 
[^-||. However, iii) has been questioned joj and ii) re- 
mained unexplained. 

In all cases where and j'^ are of comparable magni- 
tude we find that the participation number is enhanced 
by the interaction. In a mean field picture, it is the 
temporal fluctuations of the field originating from the 
partner particle which destroy the coherence necessary 
to produce localization. This explanation in terms of a 
(ij/namic-correlation-induced weakening of the influence 
of disorder can be nicely corroborated by a number of 
case studies. 

We note that contrary to previous statements, the in- 
dependence of the sign of the Coulomb interaction is 
not a general feature, but is a consequence of the im- 
posed electron-hole symmetry. In particular, displacing 
the central frequency of excitation pulses from the center 
of the absorption band, the situation changes completely. 
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Note that this choice of the excitation frequency corre- 
sponds to the reahstic situation where electron-hole pairs 
are excited close to the absorption edge in semiconduc- 
tors. Fig. 1^ shows the participation number A for light 
electrons and heavy holes, i.e., J'^ = —2J^ = 20 meV. 
The central excitation energy of the pulse is placed in the 
lower part of the pair continuum at E^xc = 40 meV. Re- 
sults averaged over 60 realizations are shown for = 80 
meV and U = 0, ±1. The results are invariant under re- 
flection of the excitation frequency through band center 
and simultaneous switching of the sign of the interaction. 
This reflects the approximate symmetry (within fluctua- 
tions in the site energy distribution) of the Hamiltonian. 

It is at first sight counterintuitive that the enhance- 
ment of the participation number is larger for attractive 
([/ = —1) than for repulsive {U — +1) interaction. This 
behavior can be attributed to the fact that for attractive 
interaction and positive masses (i.e. for excitation into 
the lower half of the excitation continuum) the electron- 
hole pair tends to stay closer together. The fluctuat- 
ing field due to the accompanying particle is then more 
pronounced as compared to the case of repulsive inter- 
action, where the mutually repulsive particle pair tends 
to be separated. Hence the dj/namzc-correlation-induced 
weakening of the influence of disorder is less effective for 
repulsive than for attractive interaction. 

Completely different behavior is found for a static fleld. 
We consider an infinitely heavy hole, j'* — 0, which now 
produces a static field, and excitation at the (interaction- 
free) band center. For both attractive and repulsive in- 
teractions, the participation number is decreased with 
respect to the noninteracting case. This result is eas- 
ily understood without invoking fluctuating fields since 
at band center electron states have maximal extent. In 
the presence of interaction, off-center states are admixed 
leading to greater confinement. The effect of the static 
interaction is thus opposite to that of a fluctuating field 

The strong retardation of the saturation in the in- 
teracting case can also be understood in our picture. 
Whether with or without interaction, the electron and 
hole wave packets spread over a range given by the single- 
particle levels involved in the optical transition just af- 
ter the short excitation pulse. The fiuctuating Coulomb 
field due to the partner particle then leads to an increase 
of the spread of the wave packets. As a consequence, 
the average fluctuating field acting on a given particle 
is reduced, which in turn tends to slow further spread- 
ing, eventually leading to the observed saturation at long 
times. The neglect of phonon interactions in our model is 
justified a posteriori. Fig. |^ makes it clear that the time 
scales between 100 fs and ~ 3 ps are fully sufficient for 
experimental observations while near-band-edge acoustic 
phonon scattering occurs on longer time scales p^ . 

Previous work on the TIP problem suggests a 

scaling of the two-particle localization length h/h ~ 



l1{U/Jf, with a = 1 or a = 0.65 |§. Our results for 
the participation number do not obey such a scaling law 
as far as the dependence on U is concerned. We obtain 
for electrons and holes, both for attractive and repulsive 
interaction, A(f/ = ±1)/A(t/ = 0) ~ K{U = Of with 
6 = 0.65 ± 0.3. For electrons and attractive interaction 
the present model predicts Ae ~ W^'^ with a larger expo- 
nent c ~ 2.2 compared to c = 1.3 for the noninteracting 
case. 

In conclusion, we have studied the localization of a 
pair of interacting particles in a situation which, in prin- 
ciple, is accessible to experiments. Optical excitation in 
the pair continuum of a disordered one-dimensional semi- 
conductor with long-range Coulomb interaction has been 
considered. Starting from a tight-binding description, 
the temporal evolution of the participation numbers of 
the electron and the hole wave packets has been calcu- 
lated by a direct solution of the equation of motion of 
the correlated material excitation within linear response 
with regard to the exciting laser field. The participa- 
tion number increases with interaction for both attrac- 
tive and repulsive interaction. We find that in general 
the degree of delocalization depends strongly on the sign 
of the interaction, in contrast to previously published 
predictions. The sign of the interaction becomes irrel- 
evant (even if the masses of electrons and holes are dif- 
ferent) only for two special situations: excitation in the 
center of the pair continuum, or excitation of the whole 
band. We have checked that this result is independent 
of the assumed form of the interaction and that it re- 
mains true also for the short-range interactions studied 
in the literature. Compared to the single-band models 
treated in the past, the present semiconductor model ad- 
mits a richer variety of phenomena, which can be qualita- 
tively explained within a mean-field picture. We empha- 
size that the enhancement of the participation number is 
clearly not due to a finite size effect, and that it should 
be experimentally observable. Ultra-short timc-of-fiight 
experiments on arrays of semiconductor quantum wires 
in the coherent limit using pump-probe techniques are a 
promising option. The enhancement should also be ob- 
servable in disordered semiconductor quantum wells. In 
this case we expect the enhancement to be even more 
pronounced, since, in contrast with one-dimensional sys- 
tems, only states close to the band edge are essentially 
affected by the disorder in two dimensions, so that the 
interaction will lead to coupling with rather extended 
states. 
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FIG. 1. Absorption spectrum of the ordered chain for 
U = 0,-1-1,-1 and for equal electron and hole masses 
{J" = —J^ — 20 meV). The spectrum of the optical pulse 
with Eexc = 80 meV and r = 100 fs is also given. 



FIG. 2. Temporal evolution of the participation number A 
after excitation by a 100 fs pulse for the interacting (solid 
lines) and noninteracting (dotted lines) cases. The disorder 
parameters are W = 20 meV and W = 60 meV 



FIG. 3. Temporal evolution of Ae (upper traces) and A^ 
(lower traces) for excitation in the lower half of the continuum 
and U = 0,±1. J" 



-23" = 20 meV, W = iO meV. 
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Brinkmann et al. : interaction assisted propagation 



Figure 1 
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Brinkmann et al. : interaction assisted propagation 
Figure 2 




Brinkmann et al. : interaction assisted propagation 
Figure 3 



